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Velocity 


There  are  multiple  plasma  models. 


o  3-Dimensions  +  3-Velocities 

o  Evolve  the  particles  position 
and  velocity 

o  e.g.  Particle-In-Cell  models 


a  Ensemble  average  of  particles 
distribution,  fs(x,\,t) 

o  Evolve  the  distribution 
function 

o  e.g.  Vlasov-Maxwell  models 
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Physical  description  of  a  fluid. 


o  Modeling  each  particle  velocity  and  position  is  not  practical. 

©  Instead  an  average  is  performed  to  give  a  statistical  description. 

o  Calculate  the  number  of  particles  per  unit  volume  having  approximately 
the  velocity  v  near  the  position  x  and  at  time  t ,  distribution  function 

/(v,x,  t) 


ps  =  ms 

psus  =  ms 


J  fs(v)d\ 
j  \fs(\)d\ 

=  P.v  =  ms  J  ww fs(\)dv,  ps  =  ^ms  J  w2fs(\)d\ 
J  WWW fs(\)d\,  hv  =  ]-ms  J  w2wfs(\)dv 


Hv  =  m 


w  =  V  -  U  v 
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Boltzmann  equation  e voles  fs, 


©  The  Boltzmann  eqn: 


dfs  dfs  qs  .  <9/v  dfs 

dt  dx  ms  [  }  d\  dt 


©  Take  the  0th ,  kst,  moments  of  the  Boltzmann  Eqn. 


v”f  dy+m>l 


n+1 


dfs  A 

dx 


+  qs  J  vn  (E  +  v  x  B)  •  ^ dx  =  ms  j  \n 


dt 


d\ 


©  Each  moment  of  the  Boltzmann  eqn  gives  an  equation  for  the  moment 
variable,  and  introduces  the  next  higher  moment  variable 

©  This  process  can  go  on  indefinitely 
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Boltzmann  equation  e voles  fs 


Ofh 

dt 


+  V  •  (ps us) 


dps 

dt 


dps  u* 
dt 


+  V  •  (/?sUsUs  +  ps I  +  11,5) 


^(E  +  UjxB)-VRJiW  + 

ms  ^ 


dps  us 
dt 


r 


des 

dt 


+  V  •  (((es  +  Ps)  I  +  n s)  •  u,  +  hs) 


—us  •  E  +  Qs  s* 

m.  z ' 


des 

dt 


r 


o  System  is  truncated  by  relating  higher  moment  variables  to  the  lower 
ones 

o  The  fluids  are  coupled  to  each  other  and  to  the  electromagnetic  fields 
through  Maxwell’s  equations  and  interaction  source  terms. 
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Advantages  of  the  Model 


LTE,  velocity  moments  ,  _____  ,  e0— M3,  me—t 0  m  ,TTTX 

Kinetic  - - - >  MFPM  — — - - — >  MHD 


Ideal  MHD  model  is  valid  when: 


o  High  collisionality,  tu/t  <C  1 
o  Small  Larmor  radius,  rLi/L  <C  1 

/  \l/2  9 

a  Low  Resistivity,  (^)  (?)  %«1 


Multi-Fluid  Plasma  Model 


o  Less  computationally  expensive  than  kinetic  models 
o  Multi-fluid  effects  become  relevant  at  small  spacial  and  temporal  scales 
o  Finite  electron  mass  and  speed-of-light  effects  are  included 
o  There  is  charge  separation  is  modeled 
o  Displacement  current  effects  are  resolved  in  the  MFPM 
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Outline 


9  Blended  Finite  Element  Method 
o  Blended  Finite  Element  Method 
o  Nodal  Continuous  Galerkin 
o  Modal  Discontinuous  Galerkin 
o  Model  Verification 


ID  Soliton  problem 

Electromagnetic  Plasma  Shock  Problem 
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The  MFPM  has  dispersive  sources. 


dQ  dV  _ 

dt  dx 


The  source  Jacobian  has  imaginary  eigenvalues 

o  Q 

The  equation  system  has  dispersive  sources 

The  dispersion  is  physical  (may  be  difficult  to  distinguish  from 
numerical  dispersion) 


9  This  dispersion  is  due  to  plasma  waves  that  result  from  ion  and  electron 
plasma  interactions  with  electromagnetic  fields 

9  An  ideal  numerical  method  for  the  MFPM  should: 

9  be  high-order  accurate  9  couple  the  flux  and  the  sources 

9  capture  shocks  o  not  impose  strict  time-step 

j|  Srinivasan  et  al,  CCP  10  (2011) 
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BFEM  SIMULTANEOUSLY  USES  CG  AND  DG. 


9  Solution  to  the  electron  and  EM  fields  is  smooth  and  does  not  shock 


9  Continuous  Galerkin 


9  Discontinuous  Galerkin 


9  Electron  fluid  and  EM  fields 


9  Multiple  ion  and  neutral  fluids 


Q  =  ^q;v« 

i 


q  =  J2cm 

i 
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Implicit  Continuous  Galerkin 


0 


o  For  this  implementation  the  balance  law  form  is  cast  as 

<9Q  dV  <9Q  _  2 

~di  +  ~dQ'fr-  S  +  kVQ“ 


®  Lagrange  polynomials  are  used  for  basis  functions,  v, 

<9Q 


b™dV=CM  >-/0 

©  # -method  time  integration 


r  dV 


+  k  /  vr\/2Qddv 
Jo 


72.(0")  =  MQ'1+ *  Q"  -  6»£r(QH+1)  -  (1  -  0)Cr(Qn)  =  0 
©  0  =  0.5  is  used  for  2nd  order  accuracy 

?(Qn)  =  ^P,  tf  (Q”)AQ  =  — 72(Q"),  Q”+1  =  Q"  +  AQ 

o|  Reddy,  An  Introduction  to  the  Finite  Element  Method  (2006) 
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Runge-Kutta  Discontinuous  Galerkin 


9Q  ,  dV  _ 

dt  dx 

o  Legendre  polynomials  are  used  for  basis  functions,  vp 
o  The  hyperbolic  equation  is  multiplied  by  the  basis  function, 

[  vp^dV  =  CJ Q)  =  [  vpSdV  -  (f  Vp^F  •  dA  +  [  •  \7vpdV 

Jn  dt  Jq  Jqq  Jq 

©  Explicit  Runge-Kutta  time  integration 
a  CFL  =  cAt/Ax  <1/ ( 2p  —  1),  p  is  the  polynomial  order 

Q*  =  Q"  +  AtCp(Qn), 

Qn+i  =  Iq*  +  Iq»  +  Ia^(q*). 

j|  Loverich  et  al,  CCP  9  (2006) 
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Convergence  of  the  BFEM. 


0 


9  Simulations  at  fixed  time- step _ 9  Simulations  at  fixed  CFL=1 
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ID  Soliton  Problem 


a  ID  soliton  is  a  two-fluid  plasma 
problem 

o  The  solution  is  smooth,  therefore 
artificial  dissipation  can  be  small 

®  The  simulation  uses  512 
second-order  elements 

o  Bz  =  1.0,  Te  =  T{  =  0.01, 

U/  =  u,  =  0 

9  ne  =  m  —  1.0  +  e-10(x-6)2 

^  Baboolal,  Math,  and  Comp.  Sim.  55  (2001) 
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DG  AND  BFEM  COMPARISON  WITH  A  SOLUTION 


—  =  1836,  —  =  1000V2,  FV  5000  cells 

csi 


®  DG  Solution  is  very  dispersive 

Hi  Hakim  et  al,  JCP  219  (2006) 


o  BFEM  is  less  dissipative  than  the 
converged  solution 
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BFEM  Computational  Cost  Savings 


0 


case 

mj/me 

c/csi 

DG  time(s) 

BFEM  time  (s) 

BFEM  cost  over  DG 

1 

25 

10/72 

0.32 

37.7 

+11681% 

2 

100 

10/72 

1.28 

37.7 

+2845% 

3 

500 

10/72 

6.82 

37.7 

+452.8% 

4 

1000 

10/72 

12.4 

38.2 

+208.1% 

5 

1836 

10/72 

23.5 

40.4 

+71.91% 

6 

3672 

10/72 

47.2 

39.2 

-16.95% 

7 

3672 

100/72 

520 

265 

-49.04% 

8 

3672 

1000/72 

5274 

2735 

-48.14% 

o  Per  time-step  explicit  DG  is  faster  than  BFEM,  but  it  requires  many  more 
time-steps 

o  BFEM  is  more  efficient  only  when  time-step  are  considerably  larger  than 
explicit  DG 
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L2-norm 


BFEM  Accuracy 


o  The  BFEM  seams  to  be  less  accurate  than  the  DG  implementation 
(~  50%) 

o  When  the  mass  ratio  is  one,  the  two  methods  have  the  same  level  of 
accuracy 

o  The  discrepancy  is  due  to  the  fact  that  the  semi-implicit  BFEM  does  not 
resolve  the  plasma  frequency  in  this  problem 
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Electromagnetic 


a  Fast  rarefaction  wave  (FR),  a 
slow  compound  wave  (SC),  a 
contact  discontinuity  (CD),  a 
slow  shock  (SS),  and  another 
fast  rarefaction  wave  (FR) 

o  The  problem  exhibits  limits 
of  MHD  and  multi-fluid 
behavior  by  changing  the 
Larmor  radius,  r l 
o  MHD:  rL  0 
o  Multi-fluid:  vL 

j|  Brio  and  Wu,  JCP  75  (1988) 
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Shock  in  Density  but  Smooth  Fields. 


o  t=0.05/cca-,  c/c„-110,  mz/mg=1836 


mass  density  Magnetic  field  (y-comp.) 

o  The  main  features  of  the  problem  are  captured  by  all  three  methods 

o  BFEM  does  not  properly  resolve  the  fast  electromagnetic  waves  which 
require  accurately  resolving  the  electron  dynamics 
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Maximum  BFEM  time-step 


A tmax  —  min 


Ax  Ax  Ax  0.1  0.1  0.1  0.1 


Ur, 


Uci  Upe  UJpi 


9  A tmax  corresponds  to  the 
maximum  value  allowed  for 
explicit  methods  based  on  the 
CFL  condition 

9  At  =  42.9A tmax  is  the 
maximum  time  step  allowed 
by  the  BFEM  due  to  ion 
dynamics 
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Effects  of  Artificial  Dissipation 


0 


9  Varying  the  artificial 
dissipation  on  the  electron 
fluid,  Ke 

9  Wave-like  behavior  of  the 
problem  is  better  resolved 

9  Amplitude  of  the  compound 
wave  increases 

9  Right  fast  rarefaction  wave  is 
not  visible 
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Effects  of  Artificial  Dissipation 


0 


o  Varying  the  artificial 

dissipation  on  the  EM-field, 

kem 

a  There  is  better  agreement 
with  the  DG  solution 

o  This  reinforces  the  point  that 
the  wave-like  behavior  arises 
from  the  interaction  of  the 
electron  fluid  with  the 
electromagnetic  fields 
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ICF  Species  Separation 


o  Deuterium  and  tritium  are 
heated  and  compressed  to 
fusion  conditions 

o  The  compression  is 
laser-driven 

o  Deuterium  can  accelerate 
faster  than  the  tritium 

o  Low  neutron  yield 

measurements  point  towards 
separation 

o  The  phenomenon  is  not 
captured  by  single-fluid 
plasma  models 
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iuo6I0l.  xu 


Accuracy 


t=150ps 


Electron  CFL=1  Electron  CFL=20 

©  The  ion  species  separation  in  both  cases  is  the  same  although  the 
solution  behind  the  shock  fronts  differ 

H  Bellei  et  al„  PoP  20  (2013) 
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Cost  as  a  fraction  of  CFL=1 


Cost  vs.  Accuracy 


o  The  solution  error  increases  as  the  computational  time  decreases. 
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L2-norm 


o  The  blended  finite  element  method  (BFEM)  is  presented 

o  DG  spatial  discretization  with  explicit  Runge-Kutta  time  integration  Ions 
and  neutrals 

o  CG  spatial  discretization  with  implicit  Crank-Nicolson  time  integration  for 
the  electrons  and  EM  fields 
o  DG  captures  shocks  and  discontinuities 
o  CG  is  efficient  and  robust  for  smooth  solutions 

o  Physics-based  decomposition  of  the  algorithm  yields  numerical  solutions 
that  resolve  the  desired  timescales 

o  DG  method  takes  less  computational  time  to  advance  the  solution  by  one 
time-step,  however  At  is  much  smaller  than  that  of  the  BFEM 

o  Computational  cost  savings  using  the  BFEM  will  only  occur  for 
relatively  large  implicit  time-steps  compared  to  explicit  time-steps 


Thank  You. 
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